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Abstract 

> 

We discuss the structure of Yukawa couplings in D-brane models in which the SM 

fermion spectrum appears at the intersections of D-branes wrapping a compact 

space. In simple toroidal realistic examples one can explicitly compute the Yukawa 
(N ■ 

couplings as a function of the geometrical data summing over world-sheet instan- 
ton contributions. A particular simple model with a N = 1 SUSY spectrum and 
' , three quark-lepton generations is studied in some detail. Remarkably, one can re- 

M-l, 



produce the observed spectrum of quark masses and mixings for particular choices 
of the compact radii and brane locations. In order to reproduce the smallness of up- 

•i— I , 

and down-quark masses branes should be located in simple geometric configurations 
leading to some accidental global symmetries. We also find that the brane configu- 
rations able to reproduce the observed data may be considered as a deformation (by 
brane translation) of a configuration with Pati-Salam gauge symmetry. The origin 
of CP- violation in this formalism is quite elegant. It appears as a consequence of 
the generic presence of U(l) Wilson line backgrounds in the compact dimensions. 
One can reproduce the observed results for the CP-violation Jarlskog invariant J as 
long as the compact radii are of order of the string scale. 

DISCLAIMER: This paper is going to be substantially revised. Althought the 
physics and general concepts are still valid, the Yukawa couplings of the particular 
model presented in this paper have a simpler form than discussed here, as we recently 
pointed out in hep-th/0302105. A properly revised version will be eventually sent as 
the paper is appropriately corrected. 



1 Introduction 



One of the most outstanding puzzles of the standard model (SM) is the structure of 
fermion masses and mixing angles. The masses of quark and leptons are clearly not 
random, showing a hierarchical structure (see table 1) with masses differing by several 
orders of magnitude. Something similar happens with the electroweak CKM mixing 



U-quarks 


u 


c 


t 




0.9-2.9 MeV 


530-680 MeV 


168-180 GeV 


D-quarks 


d 


s 


b 




1.8-5.3 MeV 


35-100 MeV 


2.8-3.0 GeV 


Leptons 


e 


M 


T 




0.51 MeV 


105.6 MeV 


1.777 GeV 



Table 1 : Masses of quarks and leptons at the Mz scale taken from the first paper in ref . [1] . 
matrix which numerically is given (for the moduli of the entries) by [2] 



\V f 



CKM\ — 



( 0.9741 - 0.9756 0.219 - 0.226 0.0025 - 0.0048 \ 

0.219 - 0.226 0.9732 - 0.9748 0.038 - 0.044 
V 0.004- 0.014 0.037- 0.044 0.9990 - 0.9993 / 



at the 90% confidence level. This is close to a unit matrix with small off-diagonal mixing 
except for the Cabibbo (12) entry which is somewhat larger. Again there is a hierarchi- 
cal structure with the third generation mixing mostly with the second generation rather 
than the first. The violation of CP is given in the "standard" parametrization [2] by the 
phase 5i3 = 59° ±13°, which may be considered rather large. A convention independent 
measure of CP-violation [3] is given by the Jarlskog invariant J which experimentally 
is [2] J = (3.0 ±0.3) x 10~ 5 . 

The understanding of the structure of fermion masses and mixings has been the 
subject of an enormous amount of effort. A phenomenological attitude widely followed 
is to consider fermion mass "textures" or definite ansatze for the form of the mass 
matrices (see, e.g., [1] and references therein). In the presence of "texture zeros" one can 
obtain relationships between the mixing angles and quark masses like the approximate 
numerical equality [4] sin 6 12 = \fn\ijm~s- Attempts to understand the presence of 
such "textures" have been made in different flavour models, one of the simplest being 
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based in Abelian flavour symmetries [5,6]. Grand unified theories combined with 
textures have also been systematically explored in order to understand the fermion 
spectrum (see, e.g., [7] and references therein). In spite of the fact that these and 
other approaches can give a semiquantitative understanding of the fermion spectrum, 
we are still lacking a full explanation of the masses and mixings in terms of a more 
fundamental theory. An obvious candidate for such a more fundamental theory is 
string theory. In fact the structure of fermion masses has been studied in a number of 
semirealistic heterotic string models. Although one can reproduce interesting features 
of the observed spectrum, one limitation of this approach is that usually the obtained 
models have (at least before field theory flat directions are taken) additional fermions, 
gauge and Higgs bosons which make the analysis complicated. 

In spite of that, string theory has a priori a theoretical advantage for addressing 
this puzzle: the Yukawa couplings may be computable as functions of the geometrical 
moduli in particular models. That is the case, for example of Abelian Zn orbifold 
compactifications in which one can use conformal field theory (CFT) techniques to 
compute them [8] . In fact, string theory contains built-in a possible mechanism in 
order to obtain hierarchical Yukawa couplings [9] . In such orbifold models some matter 
fields are localized in extra dimensions in the different fixed points of the orbifold. Thus 
one can imagine an scenario in which the Higgs field lives in one of the fixed points 
and the different right- and left-handed fermions live in different distant fixed points. 
Yukawa couplings between quarks which are distant from the Higgs are mediated by 
world-sheet instanton effects and are hence exponentially suppressed by the (distance) 2 
between the different fixed points [9] . The corresponding Yukawa couplings may be 
explicitly computed in terms of the geometry of the given orbifold. Although from the 
theoretical point of view this is fine, phenomenologically the problem arises because 
in such Ztv heterotic orbifolds searches for models with just three generations and a 
minimal Higgs sector have been unsuccessful. 

With the advent of D-branes in Type II and type I string theory, the phenomeno- 
logical possibilities of string theory have widened in several respects (see ref. [10-27] for 
the construction of explicit D-brane semirealistic string configurations). In particular, 
there is now the possibility of following a "bottom-up" approach [10,26] to the con- 
struction of semirealistic three-generation models, i.e., one can consider local D-brane 
configurations giving rise to the gauge group of the SM and three quark-lepton gen- 
erations and then embed such a local configuration in some general compactification. 
Many properties of the models depend only on the local brane configuration, rather 
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than on the details of the particular compact (e.g., Calabi-Yau) space. 

In the present paper we address the computation of Yukawa couplings in D-brane 
models. We concentrate on a recently constructed class of D-brane models [12-27] 
in which the SM chiral spectrum is obtained at intersecting branes (for short reviews 
see, e.g., [28,29]). Nevertheless, we believe the approach is much more general since 
Yukawa couplings in other type of D-brane models (like those based on D-branes at 
orbifold singularities or those involving D-branes with magnetic fluxes [30]) are related 
by T-duality and/or mirror symmetry with intersecting brane configurations [26,31]. 

One of the advantages of the present approach is that, unlike what happened 
in heterotic constructions, it is relatively easy to find D-brane string configurations 
yielding (at least for the chiral spectrum) just the particle content of the SM (or the 
MSSM) already at the string level, i.e., without any further effective field theory elabo- 
ration [16,20]. The simplest such models involve orientifold compactifications of Type 
IIA string theory on a 6-torus T 2 x T 2 x T 2 . The Yukawa couplings in these schemes 
correspond to string correlators involving the Higgs scalar, one right-handed and one 
left-handed fermion [14]. The worldsheets corresponding to those correlators have then 
a triangular shape with the fields at the intersections and the sides embedded on the 
three different intersecting stacks of branes (see figure 1). The leading contribution to 
the Yukawa couplings is then given by worldsheet instanton contributions proportional 
to exp(—S c i) = exp(—A a i> c ), where A a f, c is the area of the triangles [14]. In this pa- 
per we will focus mainly on the phenomenological applications of such computations 
and leave the derivation of general expressions for the Yukawa couplings for arbitrary 
toroidal Dp-brane configurations for a separate paper [32]. 



Figure 1: Yukawa coupling between two quarks of opposite chirality and a Higgs boson. 
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In order to show the possibilities of the approach, we compute the Yukawa couplings 
for a new simple model yielding the chiral spectrum of the MSSM at the intersection of 
D6-branes wrapping a 6-torus 1 . We provide explicit simple formulae for the relevant 
Yukawa couplings in terms of world-sheet instanton sums over the cycles of the 6-torus. 
In some particular geometrical configurations the Yukawa couplings may be written as 
products of Jacobi ^-functions with characteristics. One of the nicest features we find 
is a natural origin for the complex phases necessary for CP-violation. In the present 
approach they correspond to the generic presence of non-trivial Wilson line backgrounds 
associated to the Abelian gauge symmetries of the theory. 

We have performed a search for particular D-brane geometries in this model ca- 
pable of describing the observed quark masses, mixings and CP-violation parameters. 
We have found that for certain torus radii and D6-brane locations one can reproduce 
the observed quark masses and mixing angles. On the other hand we also find the 
interesting result that, in order to reproduce the observed size for CP- violation, the 
compact radii of two of the tori have to be of order the string scale, otherwise the 
phases get exponentially suppressed as the radii increase. The D6-brane locations re- 
producing the observed data are close to positions in the tori with certain symmetry 
properties. At those symmetric positions pairs of triangles corresponding to different 
quark generations become equal, signalling some accidental global symmetries. 

Finally we find that the observed spectrum of fermion masses and mixings is con- 
sistent with the existence of an underlying (broken) Pati-Salam symmetry [33]. More 
precisely, the brane configurations able to reproduce the observed data may be consid- 
ered as a deformation (by brane translation) of a configuration leading to a Pati-Salam 
symmetry. 

The structure of the rest of the paper is as follows. In the next chapter we give a brief 
introduction to intersecting brane models and present the specific model to be analyzed 
numerically afterwards. In chapter three we present the form of the expressions for the 
Yukawa couplings and discuss the origin of symmetries and CP-violation. In chapter 
four we perform a numerical analysis and show how all quark masses, mixings and CP- 
violation parameters may be reproduced with the given Yukawa coupling formulae. We 
also extend the discussion to charged lepton and (Dirac) neutrino masses, mixings and 
CP-violation. We leave chapter five for some final comments and conclusions. 

1 A similar analysis may be attempted with other intersecting brane models like those in refs. 
[16,20]. The present model has the advantage of having a minimal Higgs sector, J\f = 1 SUSY at all 
intersections and quite a simple geometrical configuration. This allows for rather simple expressions 
for the Yukawa couplings. 
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2 Intersecting brane standard models 



A particular class of models with intersecting D-branes has received recently consid- 
erable attention, [12, 13, 15, 16]. These are models in which one has four stacks of 
intersecting D-branes in an orientifold of Type IIA string theory. The four stacks come 
under the names of baryonic, left, right and leptonic and each stack gives rise to the 
gauge groups 577(3), SU(2)l, U(1)r and U(l)i ep t on respectively (with possibly addi- 
tional i7(l) 's). Open strings at the intersections give rise to chiral fermions with the 
quantum numbers of SM quarks and leptons. In one of the simplest class of models the 
Type IIA string is compactified on a 6-torus T 2 x T 2 x T 2 and D6-branes wrap some 
cycles of the tori. Specifically, the 7-dimensional worldvolume of D6-branes contains 
Minkowski space and the three remaining dimensions of the D6-branes wrap one-cycles 
at each of the three 2-tori of the T 6 (see figure 2). Let us denote by (n\, m l a ), i = 1,2,3 
the wrapping numbers of each brane D6 a , n l a (m l a ) being the number of times the brane 
stack a is wrapping around the rr(y)-coordinate of the i th torus. One can check that the 
number of times two branes D6 a and D6 b intersect in T 6 is given by the intersection 
number [12]: 

I ab = {n\m\ - m\n l b ){n 2 a m 2 b - m 2 a n 2 b ){n 3 a ml - m 3 a n 3 b ) (2.1) 

Open strings stretching around the intersections give rise to chiral fermions in the bi- 
fundamental representation (N a , N b ) under the gauge group of the two branes U(N a ) x 
U(N b ). Thus, these configurations yield I ab copies of the same bifundamental repre- 
sentation, providing a natural source for the observed generation replication. In fact, 
because of technical reasons 2 one performs an "orientifold" twist [34] of this theory 
by the product Q x 1Z, where Q is the worldsheet parity operator and 1Z is the reflec- 
tion with respect to the horizontal axis of the three tori. Now, for consistency of the 
construction, one has to include in the D6-brane configuration mirror branes D6 a * for 
each original £>6 a -brane, wrapping geometrical loci that are mirror with respect to the 
reflection operation 1Z. In particular, this implies that their wrapping numbers will be 
{n l a , —m l a ), rather than (n'^mj,). Being part of the same dynamical object, both D6 a 
and -D6 a * branes correspond to the same gauge group, generically U(N a ). The chiral 
fermions arising from the intersection ab*, though, will transform in the bifundamental 
representation (iV , N b ), rather than (N a , N b ). Finally, when branes a and a* coincide, 
the gauge group may be enhanced from U(N a ) to SO(2N a ) or USp(2N a ). 

2 In the example below the orientifold operation is required, e.g., in order to get a non-Abelian 
SU{2) group for left-handed weak interactions. 



5 




Figure 2: Intersecting brane world setup. Consider two D6-branes filling four non-compact 
dimensions, to be identified with M4, and wrapping three 1-cycles of T x T 2 x T 2 . In this 
example the wrapping numbers are (2, 1)(1, 0)(1, 1) (solid line) and (0, 1)(1, 2)(0, 1) (dashed 
line). The total intersection number is 2 x 2 x 1 = 4. 

Let us now show a particular D-brane configuration giving rise to three quark- 
lepton generations. This model has not been studied previously and we chose it here 
because of its simplicity and because it naturally contains a minimal Higgs sector 
which makes very simple the study of Yukawa couplings. Consider the stacks of D6- 
branes with the wrapping numbers of table 2. Genericaily, the gauge group of this 
configuration is U(3) x U(l) 3 . However, one can check that the symmetry is enhanced 
to U{3) a x SU(2)b x U(l) c x U{l)d if the brane b is located on top of its orientifold 
mirror b*. Computing the intersection numbers as above one gets the result 
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3, 


hb* 


= 3, 


lac 


-3, 


he* 


= -3 


Idb = 


3, 


hb* 


= 3, 


Idc = 


-3, 


he* 


= 3, 


he = 


-1, 


he* 


= 1, 



which corresponds to the chiral fermion spectrum of the SM (plus right-handed neutri- 
nos, see table 3). In addition there is a minimal set of Higgs multiplets if one locates 
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Ni 


{n\,m\) 


(n?,m?) 


(n?,m?) 


N a = 3 


(1,0) 


(1/P,3p) 


(l/p,-3p) 


N b = l 


(0,1) 


(1,0) 


(o,-i) 


N c = l 


(0,1) 


(o,-i) 


(1,0) 


N d = l 


(1,0) 


(1/P,3p) 


(l/p,-3p) 



Table 2: D6-brane wrapping numbers giving rise to the chiral spectrum of the MSSM. The 
discrete parameter p = 1, 1/3 describes two different sets of wrapping numbers yielding the 
same intersection numbers (2.2). The addition of the mirror branes is understood. 

the brane b on top of the brane c along the first torus (otherwise the branes do not 
intersect and the state is massive). In other words, there is a minimal Higgs sector 
with a //-parameter whose real part is given by the distance between branes b and c 
along the first torus (see figure 4 below) and its imaginary part is proportional to the 
relative phase of their Wilson lines 3 in the same torus. 
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(ab),(ab*) 
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3(3,2) 


1 
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(ac) 


Ur 


3(3,1) 


-1 


1 





-2/3 


(ac*) 


Dr 


3(3,1) 


-1 


-1 





1/3 


(db),(db*) 


L 


3(1,2) 








1 


-1/2 


(dc) 


Nr 


3(1,1) 





1 


-1 





(dc*) 


Er 


3(1,1) 





-1 


-1 


1 


(cb),(cb*) 


Hu,d 


(1,2) 





±1 





Tl/2 



Table 3: Standard model spectrum and U(l) charges. The hypercharge generator is defined 

as Qy = gQa — — ^Qd- 

If the ratios of radii in the second and third tori are equal (R^ / = R^ / R^ = 

3 By a Wilson line we mean a constant background gauge potential A a (x%) living on the worldvol- 
ume of the D6-brane a, and that depends on the three internal coordinates x" that wrap T 6 . Since 
each of these internal coordinates have the topology of a non-contractible circle, A a cannot be gauged 
away. We will take A a (x") to take values in a U(l) subgroup of the full gauge group for each stack 
a. In this way, when an open string endpoint attached to the brane a makes a closed loop in the a th 
torus, its associated wavefunction picks up a phase exp (i2i:9^). 
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X) one can check that the same M = 1 SUSY is preserved at all intersections [18,20]. So 
this configuration is (locally) M = 1 supersymmetric, and the massless chiral spectrum 
is that of the MSSM with a minimal Higgs set. In this model there are three £/(l)'s 
and only one of them (3B + L) is anomalous. As usual in this class of models (see, 
e.g., ref. [16,35]) the anomalous U(l) gets massive by combining with one RR-field. 
There are two massless £/(l)'s corresponding to (B — L) and the 3 rd component of 
right-handed weak isospin (U(l) c ). So the actual low-energy gauge group is SU(3) x 
SU(2) xU(1) B -l xU(1) c . The model may be further broken to the SM e.g., by giving a 
vev to the right-handed sneutrino. This may be triggered by the presence of a Fl-term 
for the anomalous U(l), along the lines explained in ref. [20]. We will not, however, 
discuss this issue any longer here, since our main interest will be on the computation 
of the standard Yukawa couplings of Higgs doublets to quarks and leptons. 




Figure 3: Scheme of the model in the text. Moving brane d on top of brane a one gets an 
enhanced SU (4) Pati-Salam symmetry. If in addition brane c is located on top of its mirror 
c* there is an enhanced SU(2)r symmetry. 

Note that the gauge interactions of this brane configuration may be enhanced to a 
full Pati-Salam SU(A) x SU(2) L x SU(2) R . Indeed, branes a and d are parallel (see 
table (2) )and if we put the leptonic brane on top of the baryonic branes the gauge 
group is enhanced SU (3) x U (1)b-l ~ * SU (4). Moreover if the brane c is put on top of 
its mirror c*, there is the further enhancement U(l) c — > SU(2)ji. Thus, the model we 
are studying may be considered as coming from an adjoint breaking of a Pati-Salam 
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model 4 




i = 



► 06-plane 

D6-brane a — D6-brane c 

D6-braneb •- D6-brane c* 



Figure 4: Brane configuration corresponding to the MSSM-like model described in the text, 
for the specific value p = 1. For simplicity, we have not depicted the leptonic brane nor the 
mirror a* brane. 



Before entering into the details of the Yukawa couplings, let us give a geometric 

description of the D-brane model under configuration, with the wrapping numbers given 

in table 2. Here, we will concentrate on the choice of parameter p = 1, but a equivalent 

description can be given for the model with p = 1/3. The distribution of the a, b 

and c branes is displayed in figure 4. We have not included the leptonic brane (which 

is anyway parallel to the baryonic one) for clarity. The figure shows the unit cells of 

the three 2-tori which have been chosen here square for simplicity. With an horizontal 

discontinuous arrow is denoted the "orientifold hyperpiane", i.e., the region of space 

which is left fixed under the orientifold reflection Q x 1Z with respect to the horizontal 

axis (note that the horizontal axis itself is also an orientifold hyperpiane). Note that 

the brane b, associated to electroweak SU(2)l, is taken on top of an orientifold plane so 

that branes b and b* coincide geometrically and the gauge group is enhanced from U(l) 

to USp(2) ~ SU(2). The other straight lines represent the different wrapping branes. 

Thus, e.g., the thick line represents the baryonic stack which wraps once around the 

horizontal cycle in the first torus. In the second torus it wraps three times around the 

4 In these brane configurations adjoint Higgsing corresponds to parallel brane separation and/or 
non-trivial Wilson lines along a non-Abelian subgroup. 
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vertical cycle and once on the horizontal axis. Without loss of generality we have chosen 
the cycle to pass through the origin. Finally, this baryonic stack wraps three times 
around the vertical axis (with opposite orientation) and once around the horizontal 
one in the third torus. Note that this baryonic cycle intersects three times the b, c 
and c* branes, giving rise at those intersections to the three generations of quarks. 
The left-handed quarks (ab intersections) are labelled by the index % — 0, —1, 1 and 
the origin of triplication lies in the second torus. The right-handed U-quarks (coming 
from ac intersections) are labelled by j = —1,1,0 and their triplication takes place 
in the third torus. The same phenomenon occurs for the right-handed D-quarks (now 
located at ac* intersections) which are labelled by j* = —1, 1,0. This distribution of 
intersections turns out to be important for the properties of the corresponding Yukawa 
couplings. 

Note that the branes b and c are parallel on the first torus (although they intersect 
at the other two tori). If we put them on top of each other in the first torus, massless 
chiral multiplets with the quantum number of Higgs fields of the MSSM appear. As 
we separate those two branes in the first torus, open strings have to stretch and the 
Higgs doublets get a mass rt,c/(2na'), where rj, c is the distance between those branes 
in the first torus. Alternatively, if we put branes b and c on top of each other but we 
turn on Wilson lines on the circle that they wrap on the first torus, strings between 
b and c will pick up a phase exp(2T(iipi )C ) , with ip^ = 6^ — 9^ , when going around 
such circle, so that their quantized momenta will be shifted and this will induce a mass 
ln'Phcl ' ■ Thus, both the distance r^c and the relative Wilson line (pb c are related to 
the /x-parameter in this configuration. Actually, they are proportional to the real and 
imaginary part of //, respectively. 

Notice finally that for e^ 3 ) = 0, 1/6, the brane c lies on top of its mirror c* and there 
is an enhancement U{l) c — > SU{2)r. Thus we recover a left-right symmetric model 
SU(3) x SU(2) L x SU(2) R x U(\) B -l- If in addition the leptonic brane sits on top 
of the baryonic stack (i.e., ef^ = and = e®) the gauge symmetry is further 
enhanced to a full Pati-Salam symmetry SU(A) x SU(2) L x SU(2) R . It is interesting 
to notice that in such a model with the minimal Higgs sector quark and lepton masses 
are unified and, due to the left-right symmetry the U- and D-quark mass matrices are 
proportional an hence there is no mixing. Separating the branes as in the model we are 
considering breaks the symmetry down to SU(3) x SU(2) L x U(l) 2 . Thus the distance 
of the separated branes controls how far away is the model from a PS symmetry. 
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3 The Yukawa couplings and their symmetries 



3.1 The Yukawa couplings 

In intersecting brane world configurations, chiral matter arises from brane intersections, 
being localised in a internal compact manifold. As discussed in [14], a Yukawa coupling 
arises from open strings stretching a worldsheet with triangle shape in which chiral fields 
lie at the three vertices (see also [36]). The corresponding semiclassical amplitude to 
the Effective Lagrangian has the form 

^~ exp ("w)' (3J) 

where is the area connecting the i,j and k intersections. In realistic models we 
will consider, for instance, three different stacks of branes a, b and c, such that all 
of them intersect each other, possibly several times. Let us suppose that ab intersec- 
tions yield left-handed quarks, ac intersections yield right-handed quarks and finally 
be intersections give us Higgs bosons. Then, (3.1) would give us a contribution to the 
corresponding quark- Yukawa for each triplet of intersections of each type, as depicted 
in figure 1. 

In toroidal, orbifold and orientifold compactifications, such minimal surfaces are 
given by triangles whose vertices are at the intersections and whose sides lie on the 
worldvolume of the branes. The value of such areas is then, in principle, computable 
(at least at a classical level) in terms of the geometrical data of the brane configuration. 
In fact, due to the compact toroidal geometry, for a given Yukawa coupling, there is 
more than one triangle contributing to the amplitude since there are also triangles 
wrapping a number of times the tori (see fig. 5). Thus, actually the Yukawa coupling 
comes as a worldsheet instanton sum: 

h ijk ~ £ e~ A -^^e-^M ( 3. 2) 

h,h,l3 

For a given Yukawa coupling h^k there is an infinite set of triangles with areas 
Aijk {h,h,h) with h,l2,h G Z. 5 The latter are triangles connecting copies of the 
intersections in the R 6 covering space of the 6-torus 6 . The additional phases in eq.(3.2) 

5 As we increase the values of h,l2,h, the area of the corresponding triangle also increases by a 
factor proportional to the compactification radii, and so, if these are not too small, the infinite sum 
(3.2) can be estimated by a few terms of the series. 

6 Recall that a two-torus can be defined as T 2 = R 2 /A( 2 ), where A^g = {na + mb \ n, m e Z} is 
a two-dimensional lattice that identifies points in R 2 under translations a, b. This definition can be 
generalized to e?-dimensional tori as T d = R d /A( d ), and we say that R d is the covering space of T d . 
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Figure 5: Different triangular instantons contributing to a Yukawa coupling. In a compact 
space like a T 6 , three intersection points are connected by an infinite number of triangles, 
each of them giving rise to a contribution to the same Yukawa coupling. We have illustrated 
this effect in the simple case of a two-torus (namely, the second two-torus of figure 4). In the 
figure at the right we have depicted a fundamental region of the torus, where in the left we 
have patched several copies together. There we have drawn three triangles of increasing area 
connecting the same intersection points, and thus contributing to the same Yukawa QoHqj. 

are important. As discussed in more detail in [32], 4>{hihih) are phases which are 

present in the generic case in which there are non-trivial Wilson lines of the £7(1) 's 

associated to the branes participating in the coupling They are the crucial source of 

CP-violation in the present formalism 1 . 

Let us now compute the Yukawa couplings for the specific MSSM-like configuration 

described in the previous section. For simplicity, we will suppose that branes b and c 

are one on top of each other in the first torus and that their Wilson line backgrounds 

yield equal phases in such torus. This will set the /i-term to zero. In this model the 

up-quark Yukawas correspond to triplets of intersections between branes a, b and c 

7 An additional source of complex phases may come from the inclusion of a background antisym- 
metric B-field, i.e., a complexification of the tori radii. We will not consider this possibility here. 
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(yielding abc triangles) whereas those of down-like quarks correspond to intersections 
of a, b and c* (abc* triangles). At the intersection of branes b and c, c* lie the Higgs 
fields. The left-handed quarks appear at the three intersections of a and b (labelled by 
% — 0, 1, —1) whereas the right-handed up(down)-quarks (labelled by j (j*) = 0, 1, — 1) 
live at the ac (ac*) intersections. The Yukawa couplings are thus labelled hf j: hf^, and 
depend on the areas A^-(li,l 2 ,h), Af j>f (l 1: l 2) 13). In fact, in the present case it turns 
out that all the relevant geometry takes place in the 2 nd and 3 rd torus and there is no 
instanton sum over l\. After some algebra one finds that the triangular areas are given 
by 8 [32] 

Aj (h,h) = \(^) 2 ((i? (2 V 2 ) 2 + (3 X R^r 3 ff ((i^V 3 ) 2 + (SxR^ff , (3.3) 

where we have omitted the radii subindices by defining = Recall also that 

we have x — R$ /R\^ — R^p / Ri^ ■ The r a and <f) a quantities for abc (U-quark) and 
abc*) (D-quark) Yukawas are defined in table 4. 





abc triangle 


abc* triangle 




| + e( 2 ) + l 2 




?~3 


1 + e (3) + g(3) + h 




2 


0(2) 


0(2) 


3 


0(3) + 0(3) 


0(3) _ 0(3) 



Table 4: Parameters in the MSSM-like model of table 2. 



Here i and j, j* (with values -1, 0, 1) label the intersections in the second and 

third tori, respectively, whereas l a G Z, a = 2,3, index the infinite family of instan- 

tons that contribute to the same Yukawa. The positions of the branes for quarks are 

parametrized by the continuous parameters e^ 2 \e^ and e®, their precise geometrical 

meaning being shown in figure 4. When considering leptons one has to add two more 

of these parameters and (since = e^ 3 ^) which have analogous geometrical 

meaning, exchanging branes a and d. 

8 The formula (3.3) has been computed for the specific model of table 2 with p = 1. It turns, out, 

(2) (2) 

however, that similar formulas are obtained for the case p = 1/3, under the interchange R\ <-> iij ; 

<-> . Hence, at least for the purposes of this paper, both models can be considered to give 
the same physics. 
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Concerning the phases in (3.2) one can express as [32] 

<P(k,h) = E r ^ a (3.4) 

a 

where the sum runs over the second and third tori, a = 2, 3. The phases <p a are 
linear combinations of Wilson line backgrounds of the C/(l)'s associated to the branes 
participating in the corresponding triangle. Specifically, in the present model we have 
for the U-quark Yukawa phases [32] 

2 = (-0( 2 ) - 39P) EE 9^ 

03 = (0(3) _ 30(3)) _ 0(3) = 0(3) +0(3) il5: 

For D-quarks one just replaces 9^ — > — #( 3 \ and hence 3 — + 3 ee #( 3 ) — 0( 3 ). Here 
0(«) represents a t/(l) CT Wilson line background over the a*' 1 torus (see footnote 3). 

Note that the flavour indices % for left-handed and j, j* for right-handed quarks are 
contained in r 2 and r 3 , respectively. Thus, the flavour structure depends on the brane 
configuration on the second and third torus only. Eq. (3.3) as well as the phases (3.4) 
are the basis for our phenomenological analysis in the next section. 



3.2 Symmetries and textures in the Yukawa couplings 

The Yukawa couplings so computed have certain points at which some symmetries ap- 
pear. Note to start with that the Yukawa couplings are invariant under the replacement 
(see figure 4) 

e ( 2) _^ £ (2) " n G Z (3.6) 

since this just corresponds to a relabelling of the left-handed quarks. In the same way 
one can also check the invariance of (3.3), (3.4) under 

e (3) > e (3) + k 

which is again equivalent to a relabelling of the right-handed quarks. So, in full gener- 
ality we could allow e^ 2 \e^\ to vary in between and 1/3. 

Let us now discuss certain symmetric configurations for the Yukawa couplings in 
turn. To simplify the discussion, we will define the new parameters 8^ = + e( 3 ) 
and 5^ = — e®, that appear in the expression involving the Yukawas of U-quarks 
and D-quarks, respectively. We will also take them to range in [0, 1), although that by 
the invariances (3.6) and (3.7) we can always translate any value to one in the interval 
[0, 1/3). Thus, e.g., 1/2 is equivalent to ±1/6 and so on. 
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i) ( e ( 2 \0( 2 )) = (0,0), (0,1/2), (1/2,0) 

If branes are so located, there are two identical rows in the U-quark and D-quark 
mass matrices, and hence there are one massless U-quark and one massless D-quark 
eigenvalues. This can be seen from (3.3), since for = 0,1/2 one obtains that 
Aij(h,h) = A^ij(—l 2 — 2e^ 2 \/ 3 ), and so we can exactly match both the areas and 
phases involved in the instanton sum (3.2) for i — 1 with the ones for i = —1. Thus, 
in these points there are accidental SU(2) global symmetries acting on the three left- 
handed quark multiplets. The general pattern for the Yukawa texture near one of these 
points would have the form 



h U,D = 



^ a b c ^ 

d e f 
y a b c j 



(3.8) 



where the entries are generically complex numbers. 9 It is obvious that there is, at 
least, one massless eigenstate from the presence of these two identical rows. 



ii) (5( 3 \ 0( 3 )), (#s>, 0(3)) = (0 ,o), (0,1/2), (1/2,0) 

This is somewhat analogous to the previous case. Under these conditions two 
columns of the U-quark and/or D-quark mass matrices are identical. Thus in this case 
there is an accidental 577(2) global symmetry acting on the right-handed U-quarks and 
D-quarks respectively. The mass matrix will thus have the form 



h U,D = 



a p a 
(3 a (3 
V 7 r 7 / 



(3.9) 



Again, there is at least one massless eigenstate from the presence of two identical 
columns. 



9 Strictly speaking, we find such Yukawa textures up to irrelevant phases which do not affect any 
of the physical quantities in our construction. For instance, if we consider [e^ 2 \<j>^) — (0, 1/2), then 
we find the textures 





( 


a 








\ 


( 


a 


b 


c 


\ 


h U,D = 







1 









d 


e 


f 






{ 








a 


J 


V 


a 


b 


c 


1 



which lead to exactly the same quark masses and Vqkm as (3.8) 
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iii) (e( 2 \ 0( 2 )) = (1/2, 1/2) 

In this case one can check that while A j(l 2 ,l3) = A j(— 1 2 — l,h), signalling the 
presence of two instantonic contributions of equal magnitude in the same Yukawa 
coupling, the corresponding phases associated are opposite in sign. Thus, the total 
sum of (3.2), when i — 0, can be reordered in pairs of two terms that cancel each other, 
yielding a vanishing row in both the U-quark and D-quark mass matrices 10 . In the 
same manner as in i), one can also match instanton contributions from rows i — 1 and 
% = — 1, now with a relative minus sign arising from the phase = 1/2. This implies 
the Yukawa texture 

/ a b c \ 



h u,D = 



V 



. (3.10) 

-a —b —c 

Where again a,b,c G C. Due to this very simple structure one obtains two massless 
eigenstates. 



iv) 0( 3 )), (5^, 0W) = (1/2,1/2) 

This is like the above case but case but for columns rather than rows. One thus 
gets a texture of the form 

/ a —a \ 



h U,D 



(3.11) 



p -p 

V 7 -7 J 

Again, one has two massless eigenstates. Analogous results are obtained for the leptonic 
matrices. Note that one may be close to several of this symmetry points simultaneously. 



Thus, if one is close to e.g. (S^ 3 \ (f)^) = (e^ 
texture for the D-quarks close to a pattern 

/ 



h D = 



a 







V 




-a 



{2) : 




a 



(1/2,1/2), one would have a 



(3.12) 



/ 



All these previous symmetric configurations correspond geometrically to several 

triangles getting the same area. This is translated, after summing all the instantonic 

contributions with their corresponding phases into some discrete symmetries of the 

mass matrices of the quarks, that for the above four cases may be summarized as 

10 This is particularly clear in the case with ^ = 1/3 mentioned below. In that case the Yukawa 
couplings of a row are proportional to the Jacobi theta function 8 [1/2, 1/2], which is known to vanish 
identically. 
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and whose combined effect yields one or two massless eigenvalues. Experimentally, we 
know that the up and down quarks have very small masses compared to the rest of 
the quarks. Thus in our numerical search in the next section it turns out that the 
branes will have often the tendency to sit relatively close to some of these symmetry 
points in order to reproduce the data. However, it turns out that, whenever the 
brane configuration sits simultaneously in two symmetry points, no complex phases 
arise in Vqkm and thus there is no CP violation (this can be readily seen in the case 
(S^ 3 \ (f)^) = (e^ 2 \ <fi^) = (0,0), where all the entries in (3.8) are real numbers). Note 
also that the massless modes become massive once the phases become complex, i.e., 
for 0( 2 ),0( 3 \0( 3 ) ^ 0,1/2. Let us, for instance, set e (2) = 0. Although we still have 
equal areas, that is Aij(l 2 , h) = A-ij(— l 2 — 2e^ 2 \ I3), the associated phases are now not 
equal (or opposite). This implies that the Yukawa couplings of rows % — 1 and % = — 1 
are an infinite sum of complex numbers with the same moduli but different phases. 
The interference pattern on the infinite sum of instantonic contributions in 3.2 is, thus, 
different for rows % — 1 and % — — 1, and the final Yukawas are in general also different 
complex numbers, both in modulus and phase. 



v) The point x — 1/3 

The point / = / R^ — x — 1/3 is somewhat special. Geometrically 
what happens is that branes a, d form angles ±7r/4 with the orientifold plane (branes 
b,c form angles ±7r/2). For x — 1/3 the Yukawa couplings take a particularly simple 
form. Indeed in that case one has 

1 



A(r 2 ,r 3 ) = -(2tt)' 



(R^r 2 ) +(R^r 3 



(3.13) 



Now the expression is quadratic and the whole expression for the Yukawa coupling 
(3.2) allows us to express them in terms of a product of Jacobi theta functions with 
characteristics, namely 



abc 



i/3 



,(2) 



0(2) 



(* 2 ) X 



j/3 + e( 3) + ? (3) 

0(3) + 0(3) 



(*3) 



(3.14) 



where t 2 = (^R^ /a' and t 3 = (^R^ /a', and we have imposed the data for r a , 9^ 
from table 4. The Yukawas for abc* are similar, but for the substitution (e( 3 \#( 3 )) 1— > 
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(— e®, —9^), and j h- > j*. The above theta functions are defined by 

(zy, r ) = £ e ^(«+0 a r e 2m(s+i){uH>) m (3.15) 
zez 

It is easy to show that, irrespective of the brane locations, for x — 1/3 Yukawa matrices 
all have one massive eigenvalue and two massless ones. Although the expressions for 
the Yukawas are quite elegant, in our numerical search in the next section we have not 
found parameter choices with x — 1/3 leading to a good description of the experimental 
results, although our search has not been exhaustive. 

Let us now briefly comment on some intriguing properties of the complex phases. 
One would expect that all physical quantities like mass eigenvalues and mixings should 
be invariant under an integer shift in the Wilson line phases 

Q — ► 0° + n, neZ. (3.16) 

This is not obvious from (3.4), since the coefficients r a are in general fractional. This 
implies that Yukawa matrices hfj, hf^ are not invariant under (3.16), which is quite 
puzzling. One can, however, show that this invariance is indeed present in the measur- 
able physical quantities, that is, the mass eigenvalues and the CKM mixing matrix [32] . 

One can also convince oneself of the following result concerning phases and CP 
violation. In order to get a non-negligible contribution to CP-phases, more than the 
first term should be present in the worldsheet instanton sum in eq.(3.2). If there is only 
the first term (i.e. li, l 2 , 13 — 0), then one can always reabsorb the phases coming from 
Wilson lines into quark states, so that the mass matrices are real. Note in particular 
that, if the radii are large, then the leading contribution to Yukawa couplings will come 
from the first term only, the following terms being exponentially suppressed. As a 
consequence CP-violation (i.e. the Jarlskog invariant) will be very small for large radii. 
Indeed this turns out to be the case in the numerical analysis in the next section. Thus 
we will need relatively small radii (so that several terms contribute in the instanton 
sum) in order to obtain sufficient CP-violation. 

Let us finish this section by making a remark about the absolute size of Yukawa 
couplings in this setting. The above computation corresponds to the classical contri- 
bution exp(—S c i) to the actual Yukawa correlators. In addition there are quantum 
world-sheet corrections which correspond to quantum fluctuations [8] around the flat 
triangle worldsheets. Those are expected to be flavour-independent and to affect in 
a similar way both U- and D-quark Yukawas. Thus the actual Yukawa couplings will 



•d 



5 

<f> 
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have a general form [8] Y^' = h qu h i:j ' , with h qu including the quantum fluctuation 
factor and h^- D being the classical contribution discussed above. 

Furthermore, in order to make contact with the explicit Yukawa coupling constants 
appearing in the Standard Model, all gauge and fermion fields should have kinetic 
terms normalized to one. This gives an additional factor which, e.g., in the heterotic 
case is proportional to the gauge coupling constant g. In the present case there is no 
unified gauge coupling but different couplings for each gauge group. So in our case 
one expects some geometrical average go of the gauge couplings. Finally, if we want to 
compare with fermion data at say the Z°-scale, there will be in general a loop running 
from the string scale, at which all the above expressions apply, down to the Z°-scale. 
Combining these two factors with the quantum fluctuation factor h qu one has for the 
Yukawas at the Z° scale 

Yij(Mz) = i(M z )h qu gvh tj = h h tj (3.17) 

where i{M z ) is a renormalization group factor. Those corrections are expected to be 
similar for both U- and D-quarks, since the leading effects should come from QCD 
loops. Thus we will take = = h . On the other hand, the running effect on 
leptons should be much smaller and one expects approximately to have = h /^(M z )- 
In any event, in the numerical computation we will leave h and as free parameters. 



4 Reproducing the experimental results 

Let us now see if the above brane configuration is able to reproduce the quark spectrum, 
mixing angles and CP violation. We will also consider the charged lepton and Dirac 
neutrino masses below. The Yukawa couplings we computed in the previous section 
applied at the string scale. If the string scale is high the running effects for the quarks 
down to the weak scale may be large. As we mentioned above we will reabsorb those 
effects in the definition of Hq ,d = h and h$ and then we will compare with data at 
the Z° scale. 

As we said, this model can also be understood as a Pati-Salam model in which brane 
separation (adjoint Higgsing) breaks it to SU(3) x SU(2) L x U(l) 2 . In the presence of 
exact Pati-Salam symmetry and a minimal Higgs sector one has (before r.g. running) 

M v = M v = tg(3M D =tg^M L (4.1) 

where tg(3 = ^f^- This leads to identical hierarchies for all fermions and no CKM, so 
this cannot give, as it stands, a good description of the SM fermion spectrum. However 
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this may be a good starting point since after all the experimental hierarchical structure 
is somewhat similar for the different quarks and leptons and the CKM mixing matrix 
is not far away from unity. So it seems it could be a good idea to perturb around such 
a D-brane configuration. We now separate from this PS symmetry point by setting 



e\ 2 > ^ e (2) and e} 3) + e (3) . Then the SU{4) breaks to SU(3) x (7(1). In addition we 
set e (3) ^ 0, 1/6 and then SU(2) R breaks to U(l) c . Naively, our field theory experience 
would tell us that this is not going to affect the electroweak Yukawa couplings, those 
are different couplings which have nothing to do with the adjoints. However, the 
implementation in terms of wrapping D-branes tells us that the Yukawa couplings 
are affected by this breaking: since the relative brane locations vary, the areas of the 
triangles also vary and the equations (4.1) no longer hold. 11 

Once we have separated from the PS symmetry point, the mass matrices for U- and 
D-quarks are no longer proportional and are given by 

= h hfj<H> ; M? = cotg{f3)h hf j <H> (4.2) 

respectively, where 

V2M W 



^J\<H>\ 2 + \<H>\ 2 = 



174GeV 



9l 



(4.3) 



and hf j: hf- are taken from eqs. (3.2), (3.3), (3.4) and the data on table 4. They will 
explicitly depend on the geometry of the tori (R^ 2 \ R^ 3 \ x), the e's which parametrize 
the brane positions and the phases 9 coming from the Wilson lines. The mass matrices 
are in general not symmetric and becomes convenient to work with the hermitic matri- 
ces M u (M u y and M D (M D y in order to obtain the mass eigenvalues and the unitary 
matrices U^, 



(uvy(M u (M u y)uv = 




x 



(c/£ y (M D {M D y) u 



( 



v 










2 

m c 



m\ 




\ 


m\ ) 
\ 




m\ ) 



(4.4) 



11 From the field theory point of view this means that the non-renormalizable operators of type 
if) Rip lH (adjoint)" are playing a role. The difference is that in our D-brane realization all these non- 
renormalizable terms are automatically summed to all orders by considering how the area of the 
triangles involved in the Yukawa computation have changed. Note that the e's and </>'s parametrizing 
brane locations and Wilson lines correspond to vev's for complex adjoint scalar fields in these brane 
settings. Thus, e.g., an expansion on small separations e's from the PS limit correspond to an expansion 
involving powers of adjoint scalars. 
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which diagonalize them. Then the CKM mixing matrix will be given by 

Vckm = (U^ UE (4.5) 

and will generically be complex. In the present case it is convenient to express CP- 
violation in a convention independent way, i.e., in terms of the Jarlskog invariant J [3] 
which may be computed in terms of four complex entries of Vckm, for example: 

J = Im(V cs V: s V ud V: d ) . (4.6) 

In terms of the 'standard parametrization' [2], this may be expressed in terms of three 
real mixing angles #y and the complex phase 5± 3 : 

J = s 12 s 13 S23C 12 C23C 2 13 sin(5 13 ) (4.7) 

where Sij = sm(%), = cos(%). 

The mass matrices depend on the following 8 real parameters: h , tg[3, F6 2 \ x, 
e^ 2 \ and e^. In addition there are the phases 8^ , 9^ and 9^ 3 \ We have made a 
computer search for choices of these parameters able to describe the observed spectrum 
of quark masses, mixing angles, CP-violation phase and charged lepton masses. We 
will discuss first the quark sector. We find different regions in parameter space in which 
all those data can be reasonably adjusted and we show some possible choices in table 5 
and the corresponding results in table 6. Note that the values for ho and tgf3 are fixed 
by the condition that they reproduce the results m,b(M z o) = 3.0 GeV, m t (M z o) = 180 
GeV. 



Model 


ho 


tgp 






X 


e (2) 


JO) 


J(3) 




0( 3 ) 


$3) 


A 


1.169 


33.33 


0.105 


0.66 


2.552 


0.09 


0.50056 


0.03 


0.43 


0.43 


0.485 


B 


0.715 


74.48 


0.245 


0.17 


1.71 


0.024 


0.5004 


0.041 


0.1405 


0.472 


0.384 


C 


0.715 


74.20 


0.2516 


0.1714 


1.71 


0.024 


0.50059 


0.029 


0.14 


0.47383 


0.386 


D 


0.493 


11.04 


0.1 


0.15 


1.75 


0.452 


0.0011 


0.5091 


0.4539 


0.0 


0.25 


E 


0.597 


63.7 


0.3024 


0.2 


1.75 


0.186 


0.0009 


0.1254 


0.098 


0.50001 


0.403 



Table 5: Some choices of parameters leading to the results shown in table 6. The radii 
are given in units of M~ x . 

Let us summarize our findings on a few points: 
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Model 


m u 


m d 


m s 


m c 


m h 


m t 


|Vi 2 | 


1^23 1 


1^13 1 


J 


A 


1.28 


4.32 


81.75 


669.1 


3.0 


180.0 


0.224 


0.0499 


0.0023 


0.9 x 10~ 5 


B 


4.37 


6.90 


65.73 


576.46 


3.0 


180.0 


0.2194 


0.0386 


0.0049 


3.0 x 10~ 5 


C 


2.81 


7.90 


71.07 


646.97 


3.0 


180.0 


0.2259 


0.0408 


0.0048 


2.88 x 10~ 5 


D 


0.95 


1.97 


73.03 


607.88 


3.0 


180.0 


0.2199 


0.0484 


0.0025 


2.61 x 10~ 5 


E 


0.91 


1.04 


104.05 


709.57 


3.0 


180.0 


0.225 


0.0381 


0.0058 


1.14 x 10~ 5 



Table 6: Quark masses, mixing angles and Jarlskog invariant at the Z° scale for the 
model parameters of table 5. The masses for the two heaviest quarks are given in GeV 
and the rest in MeV. 

• Appropriate results are obtained if the branes are close to some symmetry points 
as discussed in the previous section. We have found all the acceptable results to 
be very close to points with 

5( 3 ) = 0( 3 ) = 0, 1/2 (4.8) 

with a good accuracy. Also some examples are in addition close to other symme- 
try points like e^ 2 -* = (p^ = 1/2 (example D) . It is easy to understand the origin 
of eq.(4.8). The biggest intergeneration hierarchy in the standard model is that 
of m up /mtop oc 10~ 5 (compared to the D-quark hierarchy which is md/mb oc 10~ 3 
only). If the brane b is close to a symmetry point 5^ = 0, 1/2 in the third torus, 
then we mentioned in the previous section that a massless U-quark eigenvalue is 
obtained, explaining the presence of the U-quark hierarchies. On the other hand, 
the hierarchy of the D-quarks is milder and is obtained by having, e.g., either <r 3 ) 
or/and ~ 0, 1/2. 

It is interesting to show the type of texture for the U- and D-quark masses that 
one obtains. For example for the choice of parameters D in table 5 one finds 
for the moduli of the entries \Y^' D \ and the phases a^- D (in radians) of Yukawa 
couplings 





I 0.958 


0.843 


0.959 N 




( 


-0.290 


-0.288 


-0.290 N 




1^1 = 


0.245 


0.221 


0.246 






-0.455 


-0.453 


-0.455 


■ (4-9) 




v 0.783 


0.697 


0.784 j 




V 


0.504 


0.502 


0.504 j 






/ 0.185 


0.087 


0.195 \ 




f 


-0.494 


-0.252 


-0.086 \ 




\Y D \ = 


0.056 


0.032 


0.058 






-0.614 


-0.454 


-0.282 


. (4.10) 




v 0.162 


0.086 


0.170 , 






0.300 


0.509 


0.669 j 
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8 = 1/2 8 = 1/4 




i=0 i=-l i = l 

► 06-plane 

» D6-brane a D6-brane c 

► D6-brane b D6-brane c 



Figure 6: Brane configuration corresponding to the symmetry point in (4.11), (4.12) and 
(4.13). For concreteness, we only show the 2 nd and 3 rd tori. 

One can qualitatively understand this structure by noting that the brane location 
parameters D are relatively close to a symmetry point with 

e ( 2 ) = 0(2) = 1/2 (4.11) 

6 (3) = ^(a) = o ( 412 ) 

<5 (3) = 1/2 (4.13) 

and <fr^ at some intermediate value between and 1/2 (namely ( - 3 - ) = 1/4). 
The brane configuration in the 2 nd and 3 rd tori corresponding to this symmetric 
configuration is depicted in fig. 6. According to our discussion on symmetries 
and textures in previous section, due to the proximity to the condition (4.11) 
both textures of U-quark and D-quark Yukawas should present a behaviour quite 
similar to (3.10), with the entries of one row much smaller than the other two. 
In addition, due to (4.12), the U-quark matrix should have two equal rows, as 
in (3.9), whereas the D-quark matrix should have some intermediate behaviour 
between (3.9) and (3.11). Indeed, the mass matrices in (4.10) behave qualitatively 
like that. So the example D in table 5 is able to reproduce the experimental results 
due to its proximity to points with special symmetries. Something analogous 
happens with the other choices of parameters in table 5, being all relatively close 
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to some symmetry point 



• The brane configurations able to reproduce the data may be understood as a 
continuous deformation of a left-right symmetric configuration. As we said, for 
g( 3 ) = (or, equivalently, 5^ = 5^) the brane c is on top of its mirror c* and 
there is a gauge enhancing U (l) c — > SU (2)# to a left-right symmetric model. Un- 
der those circumstances the U-quark and D-quark mass matrices are proportional 
and hence there is no CKM mixing. As we switch on e^ 3 - 1 7^ Mu and Md cease 
to be proportional and CKM mixing is generated 13 . This may be numerically 
seen in fig. 7. 

In fig. 7(d), the solid and dashed lines show the masses of U- and D-quarks masses 
respectively of the three generations as reproduced by model D in table 5. One 
can restore SU(2) R symmetry by setting 5^ = 5^ (fig. 7(c)) and setting 0® = 
0® (fig. 7(b), (a)). Then one observes that the masses of U-quarks and D-quarks 
are proportional, as it should in a left-right symmetric model. One can also 
check that in this LR-symmetric limit the mixing disappears. Note also that one 
observes that this SU (2) R breaking is the cause for the observed inverted hierarchy 
m up < rridown as well as the appearance of mixing (and also CP- violation). 

• One can easily obtain fermion mass hierarchies consistent with experimental 
data. However, getting substantial CP violation requires that there are several 
instanton numbers contributing significantly to the world-sheet instanton sum 
in eq.(3.2). This happens when the torus radii and are smaller than 
M" 1 , 14 otherwise only the first term contributes and CP- violation is suppressed. 
Note however that the radii R^ and R^ cannot be too small, since then the 
instanton sum diverges and the quantum piece of the Yukawa coupling becomes 
relevant. We have checked in all the numerical examples that the corresponding 
instanton sum converges. In most of the cases the terms in the sum eq. (3.2) 
become negligible for |/ 2 |, \h\ > 6. 

12 They cannot, however, sit exactly on such points, since this would imply no CP violation. 

13 One can measure how strongly SU (2) r is broken in terms of 2e^ — 5^— d^ 3 \ which is the distance 
between branes c and c* in the third torus as well as the Wilson line background 9^ = 8^ — 9^ 
breaking SU(2)r. In fact the mass of the Wr gauge boson may be written in terms of those two 
parameters. Going through the examples of the tables one finds that Mw R is always of order the 
string scale. So one can conclude that the underlying SU (2)r symmetry is strongly broken. 

14 Note that this does not affect the radius of the first torus whose size is unconstrained by Yukawa 
couplings. 
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Figure 7: The masses of U-quarks (solid line), D-quarks (dashed line) and charged leptons 
(dashed-dotted line) at the Z° scale for the three generations. Figure (d) shows the results 
for the choice of parameters D in tables 5 and 7, in good agreement with experiment. Figures 
(a), (b) and (c) show how that brane configuration may be continuously reached from an 
underlying Pati-Salam configuration by shifting of branes and Wilson line breaking (see text). 

• The complex phases coming from Wilson lines not only may account for the 
observed CP-violation. In addition, they have important influence on the values 
of the masses of the lightest fermions. Recall that, e.g., as we mentioned in the 
previous section the symmetry points with e^ 2 \ S^ 3 \ 6^ = 0, 1/2 give rise 
to massless fermions only if, in addition, one has real phases <jr 2 \ (p^, <fi^ = 
0, 1/2. Thus, for non-vanishing complex phases these massless modes will get 
massive. 
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We conclude that one can reproduce in a rather satisfactory way the quark masses, 
mixings and CP-violation in terms of the present model. Note that we have not done 
a fully systematic search for ranges of parameters able to reproduce the experimental 
results and the choices shown in the tables are just possible examples. Some values 
of parameters are particularly successful, like the choice D in table 5. The values we 
found for tgf3 vary in the range 10 — 74 and the value for ho is of order one, as expected 
from our discussion at the end of section 3. 

Although there are a number of free parameters, note that the brane location pa- 
rameters e, 5 and phases in almost all cases sit close to symmetry points. One may 
speculate that those parameters could be fixed by some symmetry in the dynamics 
which eventually determines the geometry of the brane configuration. One can con- 
sider an scenario in which all the branes are located on top of some symmetry point 
like, e.g., that in (4.11), (4.12), (4.13) and some (e.g., loop) corrections give rise to 
a small deviation from those values. In such scenario the number of free parameters 
would be essentially reduced to the geometric moduli and tgfi. 

A similar study can be done for charged leptons and Dirac neutrino masses. In fact, 
once fixed and 5^ (plus phases) to reproduce the charged lepton spectrum, Dirac 
neutrino masses as well as Dirac mixing angles and leptonic CP-violating phases are 
fixed, since the geometry of the whole brane configuration gets completely determined. 
This is a very interesting point. However, in the case of neutrino masses the structure 
of the Majorana mass for the right-handed neutrinos is very relevant. As we mentioned 
above, once some linear combinations of sneutrinos get a vev (which may be triggered 
by a Fl-term), the SU(3) x SU(2) x U(l) c x U(1) B -l symmetry is broken to the SM 
group. One can also see that in the same process r.h. neutrinos get a mass [20]. 
The latter mechanism is somewhat different to the class of Yukawa couplings we are 
considering here. Thus we will present a complete analysis of the leptonic and neutrino 
sector in this class of models in a separate publication. Here we will contempt ourselves 
with the analysis of the charged leptonic sector. Again one can reproduce the observed 
spectrum of charged leptons in terms of the location parameters of the leptonic brane 
d, ef^ and (and phases). In addition, once we have adjusted the quark and charged 
lepton sector, the Dirac neutrino masses and (Dirac) neutrino mixing matrix as well as 
leptonic CP-phases are fixed. Table 7 shows some choices of the leptonic parameters 
(corresponding to the previous quark sector choices A-E) giving rise to the results in 
table 8. Of course, only the charged lepton masses are directly measured at present, 
and those can be adjusted remarkably well. As we said, the rest of the information 



26 



is certainly relevant for the structure of physical neutrino masses and leptonic CP 
phases, but a knowledge of the r.h. Majorana masses should be first provided. Let 
us only mention that the leptonic CP-phases may be large in some cases and also 
that sometimes one of the Dirac neutrino masses may be essentially vanishing. That 
happens, for instance, in model D because one sits on the point (#} 3 \ 0f ) = (1/2, 1/2), 
giving a vanishing column in the Dirac neutrino mass matrix. The corresponding 
sneutrino may then acquire a large vev (breaking (B-L)) without giving rise to a large 
mass term of type HL. This large vev would then give rise to large Majorana masses 
for all r.h. neutrinos, and a standard seesaw mechanism could then be at work. 

We also find that the brane configurations able to reproduce in addition the observed 
spectrum of charged leptons may be understood as a deformation of a Pati-Salam 
symmetric configuration. As we said, if the leptonic brane sits on top of the baryonic 
stack (i.e., ef^ = and = e®) the gauge symmetry is further enhanced to a full 
Pati-Salam symmetry SU(A) x SU(2) L x SU(2) R . Numerically we find (modulo 1/3) 
e ( 2 ) _ e ( 2 )~ ^( 3 ) _ §(3) ^ o.04 — 0.09 indicating that the brane configuration able to 
reproduce the quark masses and mixings as well as the charged lepton masses may be 
considered as a deformation of a brane configuration leading to a Pati-Salam SU(4) 
symmetry 15 . This is again illustrated in fig. 7. In fig. 7(a) we have set e*- 2 -* = ef\ = 

(3) 

e) and also equal phases. Then the masses for D-quarks and charged leptons become 
proportional. In fact in the PS limit they should be equal at the string scale, but recall 
that we have allowed for different original values for ho and h l , which could take in to 
account e.g., the renormalization group running from the string scale to the Z° scale. 



Model 


hi 


4 2) 




J. 2 ) 


A 3 ) 


A 


0.684 


0.35 


0.0011 


0.35 


0.3104 


B 


0.43 


0.4 


0.155 


0.23 


0.29 


C 


1.51 


0.17 


0.03 


0.47 


0.38 


D 


0.166 


0.499 


0.508 


0.1845 


0.25 


E 


0.141 


0.0 


0.06 


0.06 


0.14 



Table 7: Choices for the leptonic parameters corresponding to the ones in table 5 and 
leading to the results in table 8. 

15 Note that this applies if the baryonic and leptonic branes sit on top of each other in the first torus, 
which is not necesarily the case. The Yukawa couplings are not sensitive to the geometry in the first 
torus. 
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Model 


m e 




m T 








\vh\ 


\v&\ 


1^3 1 


Ji 


A 


0.512 


105.66 


1777 


18.15 


3.615 


124.77 


0.127 


0.098 


0.0002 


1.16 x 10~ 6 


B 


0.516 


105.59 


1777 


309.5 


4.683 


83.95 


0.052 


0.228 


0.021 


15.3 x 10~ 5 


C 


0.513 


105.48 


1777 


1.90 


1.698 


160.22 


0.021 


0.057 


0.0018 


5.7 x 10~ 7 


D 


0.516 


105.52 


1777 


0.0 


0.652 


229.45 


0.200 


0.0648 


0.0232 


5.9 x 10~ 7 


E 


0.511 


105.56 


1777 


35.72 


5.989 


82.27 


0.009 


0.128 


0.0044 


1.39 x 10~ 6 



Table 8: Charged lepton and Dirac neutrino masses, mixing angles and (leptonic) 
Jarlskog invariant at the Z° scale for the model parameters of tables 5 and 7. All 
masses are given in MeV except for the two heaviest (Dirac) neutrino masses which 
are given in GeV. 

5 Final comments and conclusions 

A number of comments are in order: 

• A few words about the particular brane configuration analyzed. For the compar- 
ison with experimental data we have considered a new D6-brane configuration 
which is particularly simple. At the local level the model has Af = 1 SUSY and 
the MSSM chiral spectrum. It has the minimal Higgs sector of the MSSM and 
the /^-parameter is governed by both the distance between the branes b and c and 
their relative Wilson line in the first torus. The gauge group is that of the SM 
plus an extra U(l) familiar from left-right symmetric models. The latter may be 
broken spontaneously by a right-handed sneutrino vev, as explained in ref. [20]. 
If this configuration of D6-branes is wrapping a 6-torus, there are in addition 
adjoint scalars and fermions with respect to the gauge group, which would be- 
come massive after SUSY breaking. Let us remark that Type II string theory 
models with D-branes like this have consistency constraints from the cancella- 
tion of the Ramond-Ramond (RR) tadpoles. Specifically, the theory has certain 
antisymmetric fields under which D-branes (and orientifold planes) are charged. 
Tadpole conditions come from imposing that the overall RR-charges of the D- 
brane configuration vanishes. In the example considered in the text the D-brane 
configuration does not cancel by itself the RR-tadpoles. Thus in order to cancel 
RR-tadpoles there should be further D-branes in the model. Since however, the 
SM brane configuration is anomaly-free, the extra branes will have vanishing in- 
tersection numbers with the SM branes, i.e., they will not give rise to extra chiral 
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massless modes, being some sort of 'hidden sector' in the model [20]. One can 
also check that in the simple toroidal case the complete D-brane configuration 
including the SM branes and the extra D-branes is necessarily non-SUSY [12,20]. 
Thus one would be forced to lower the string scale (by making some compact 
volume large [39] ) in order to avoid a gauge hierarchy [13, 26] . Alternatively 
one could consider fully Af = 1 SUSY extensions of this model, which is any- 
way Af = 1 supersymmetric locally. In particular it would be interesting to see 
whether one can embed this very simple brane configuration into a fully Af = 1 
SUSY model as in refs. [18,27]. 

• Irrespective of its construction as a string model, one may consider the brane 
configuration itself, with the fermions and Higgs fields at intersections, as a phe- 
nomenologically viable model for the understanding of quark and lepton masses 
and CP violation. The idea would be that analogously simple D-brane configura- 
tions incorporating similar mechanisms and symmetries could possibly be present 
in large classes of models with D6-branes wrapping cycles on compact manifolds 
(like, e.g., a Calabi-Yau). Note in this respect that the size of the string scale, 
either high or low, has no relevance for the computation of the Yukawa couplings. 
Thus our Yukawa coupling formulae may be applicable to both low string models 
and more conventional models with the string scale of order the GUT scale. We 
have seen, however, that in this toroidal example the size of the second and third 
tori have to be of order of the string scale in order to adjust the fermion spectrum 
data. 

• The explicit formulae for Yukawa couplings provide a new laboratory for the un- 
derstanding of the structure of fermion masses. Some properties of these formulae 
include possibilities phenomenologically explored previously. Thus for example, 
our Yukawa couplings depend on the brane location parameters, the e's. Those 
correspond to vacuum expectation values of adjoints under the full gauge group 
(including C/(l)'s). Expanding our formulae on the e-parameters would give rise 
to non-renormalizable couplings of the form Hip R ip L (adjoint) n . This kind of cou- 
plings have been used since a long time in attempts to understand the SM Yukawa 
couplings. Unification symmetries like GUT's have also played an important role 
in previous analysis. In the specific example studied in detail a Pati-Salam sym- 
metry SU(A) x SU(2) L x SU(2) R is underlying the brane configuration. The 
present scheme has also some analogies with the structure of Yukawa couplings 
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in heterotic Z^v orbifold models [8,9], in which fermions and Higgs multiplets lo- 
cated at different fixed points interact via world-sheet instanton effects 16 . This is 
also somewhat analogous to recent attempts [37] to understand the fermion spec- 
trum in terms of extra dimension models in which the fermions have different 
locations (with, e.g., a Gaussian wave function) and the different wave function 
overlap with the Higgs field determines the structure of Yukawa couplings. 

• Other property that we find is the presence of accidental (approximate) global 
(e.g., SU(2)) symmetries in the structure of Yukawa couplings 17 . For some geo- 
metrically symmetric brane locations some of the triangles formed by the branes 
become identical. At those locations (and for real phases) two columns and/or 
rows of the fermion mass matrices become identical, signalling the presence of 
a massless fermion. The numerical analysis shows that in order to have suffi- 
ciently light up- and down-quarks the branes should be close to these symmetry 
points. We have also found that for other symmetric brane configurations there 
are cancellations between different instanton contributions giving rise to vanish- 
ing columns and/or rows in the mass matrices. It is the interplay of these two 
effects which turns out to be able to reproduce the observed fermion spectrum, 
rather than any exponential suppression of the lightest fermions. 

• One of the nice insights of the present construction is that of the origin of CP- 
violation. In our constructions complex phases appear associated to the generic 
presence of non- vanishing Wilson lines of the £7(1) 's of the model. We believe 
this is rather general and is not specific to our particular constructions. This 
means that the origin of CKM CP- violation is an open string effect, i.e., it is 
associated to the D-branes, rather than to the bulk fields. This is interesting 
because if promoted to a fully Af = 1 setting, it may provide a solution to the 
supersymmetry CP-problem, that is, why the phases of soft terms are so small 
(which is required from EDMN limits) compared to the CKM phase (for some 
recent work on CP- violation in string models see, e.g., [40] and for the same issue 
in extra dimensions see [41]). If soft terms physics comes from the bulk (closed 
string) gravitational sector, they may be all real without affecting CKM mixing. 

16 From the string theory point of view there are a number of differences though. For example, in 
the present case the Yukawa couplings come from open strings in the disk. In the heterotic string 
there are only closed strings and the relevant amplitude appears on the sphere. 

17 The presence of flavour 5/7(2) global symmetries in the mass matrices have been previously 
considered, e.g., in [38]. 
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In other words CKM CP-violation would come from the open string sector and 
soft-term CP-violation from the closed string (bulk) sector. 

In summary, we have shown how simple configurations of D-branes wrapping a 
compact space may give a good quantitative description of quark masses and mixing 
angles as well as CP-violation and charged lepton masses. We believe that the approach 
is more general than the particular example numerically explored since all realistic D- 
brane models constructed up to now may be re-expressed as coming from D6-branes 
wrapping cycles in some compact space [26]. 

One of the nice features of this approach is that one can obtain simple explicit 
formulae for each Yukawa as a sum over world-sheet instanton contributions. A non- 
vanishing value for U(l) Wilson lines lies at the origin of complex phases and CP 
violation. Symmetric configurations of the branes lead to approximate symmetries of 
the fermion mass matrices, giving rise to particular patterns in the fermion spectrum. 
Interestingly enough our analysis shows that the experimental data is consistent with 
the existence of an underlying Pati-Salam symmetry which is broken by shifting of the 
brane locations (adjoint Higgsing in the field theory language). This shifting is also the 
cause of the reversion of the up-down hierarchy (i.e., m up < m down ), the presence of 
CKM mixing and CP violation. It is certainly worthwhile to study the neutrino mass 
spectrum from this perspective. We hope to report on this issue in the near future. 
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